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Abstract 

In Universal Extra Dimension models, the lightest Kaluza-Klein (KK) particle is generi- 
cally the first KK excitation of the photon and can be stable, serving as particle dark matter. 
We calculate the thermal relic abundance of the KK photon for a general mass spectrum 
of KK excitations including full coannihilation effects with all (level one) KK excitations. 
We find that including coannihilation can significantly change the relic abundance when the 
coannihilating particles are within about 20% of the mass of the KK photon. Matching the 
relic abundance with cosmological data, we find the mass range of the KK photon is much 
wider than previously found, up to about 2 TeV if the masses of the strongly interacting level 
one KK particles are within five percent of the mass of the KK photon. We also find cases 
where several coannihilation channels compete (constructively and destructively) with one 
another. The lower bound on the KK photon mass, about 540 GeV when just right-handed 
KK leptons coannihilate with the KK photon, relaxes upward by several hundred GeV when 
coannihilation with electroweak KK gauge bosons of the same mass is included. 



1 Introduction 



One of the most important astrophysical challenges is to understand the nature and identity 
of dark matter. Perhaps the most interesting candidate is a neutral, stable, weakly interacting 
particle arising from physics beyond the Standard Model. This is consistent with a wide range 
of data, such as measurements of the cosmic microwave background (CMB) radiation that imply 
the abundance of dark matter is roughly a factor of 6 times the abundance of baryons as 
well as comparisons between large scale structure simulations and galaxy survey observations 
that suggest the bulk of dark matter is cold [|| . 

Universal Extra Dimension (UED) models || (see also @), in which all of the Standard 
Model fields propagate in extra dimensions, provide an interesting example of Kaluza-Klein 
(KK) dark matter §,§,§• This is because bulk interactions do not violate higher dimensional 
momentum conservation (KK number), and in these models all of the couplings among the 
Standard Model particles arise from bulk interactions. To generate 4D chiral fermions, the 
extra compact dimension(s) must be modded out by an orbifold. For five dimensions this 
is S /Z%, while in six dimensions T 2 /Z2 is suitable and has other interesting properties |^] 
including motivation for three generations from anomaly cancellation Q and the prevention of 
fast proton decay Q. An orbifold does, however, lead to some of the less appealing aspects of 
the model. Brane-localized terms can be added to both orbifold fixed points that violate KK 
number. If these brane localized terms are symmetric under the exchange of the two orbifold 
fixed points, then a remnant of KK number conservation remains, called KK parity. All odd- 
level KK modes are charged under this discrete symmetry thereby ensuring that the lightest 
level-one KK particle (LKP) does not decay. The stability of the LKP suggests it could well be 
an interesting dark matter candidate. 

We calculate the thermal relic abundance of dark matter including coannihilation with all 
level-one KK particles. This is an important calculation since there are many level-one KK 
particles that are near enough in mass to the LKP that they are expected to be important 
to determine the relic abundance. The identity of the LKP depends on the mass spectrum of 
the first KK level. We assume it is the lightest KK photon. This has been widely advocated 
for both its properties as dark matter and is favored by the contributions to the masses of the 
level one KK particles from one-loop radiative corrections || [7J. The spectrum of level-one KK 
particle masses is also not known. However, certain finite and log-enhanced one-loop radiative 
corrections |l(J to the masses have been calculated in Q. The size of these corrections depend 
on two unknowns: the size of the log, i.e., the cutoff scale of the theory, and the size of the 
matching corrections evaluated at the cutoff scale. Generically one would expect the calculable 
log-enhanced corrections to dominate over the uncalculable matching corrections. However, the 
log is not particularly large relative to the finite matching corrections. Indeed, should the cutoff 
scale be close to the scale of the KK particle masses, uncalculable matching corrections would 
be of order the calculable log-enhanced corrections, and then one would like to know the relic 
abundance for a general spectrum. This is the calculation we carry out in this paper. 

We present our results with three scenarios of level-one KK particle mass spectra. First, 
we consider coannihilation with all level-one KK particles that carry electroweak, but not color 
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charges. Next, we consider coannihilation including colored particles. In both of these cases, we 
take all of the level one KK particles (other than the LKP) to have the same mass m (1) (1 + 5). 
Finally, we consider coannihilation with all particles with a mass spectrum derived from the 
loop corrections from ||. We believe these three scenarios provide a good cross section of the 
effects of coannihilation. We also provide formulae for the (co) annihilation cross sections so that 
any other more complicated spectra could easily be calculated based on our results. 

The format of this paper is as follows. In Sec. Q we discuss KK dark matter including past 
results and motivate the need for the calculations we carry out. In Sec. |3|, we briefly describe 
the method used to calculate the relic abundance of the KK photon, and the effect of including 
coannihilating particles in this calculation. In Sec. we consider the effect of coannihilation on 
the relic abundance, and use this to obtain limits on the mass of the corresponding dark matter 
candidate, the KK photon. We find mass ranges in three scenarios: (1) when electroweak 
KK excitations coannihilate, (2) when all KK excitations coannihilate, taking all KK particles 
(except for the LKP) to have the same mass, and (3) when all KK excitations coanihilate 
assuming the mass spectrum is determined by the one-loop radiative corrections in j|] with zero 
matching corrections. Appendices |X| and ^ reviews the Lagrangian and Feynman rules in UED. 
The diagrams and cross sections pertinent to these results are summarized in Appendices |C] and 

0- 



2 KK dark matter 



The relic density of the KK photon was first calculated by Servant and Tait ||. Using their 
results, one finds that its thermal relic abundance of the KK photon matches the WMAP 
observations for cold dark matter when the KK photon mass is between about 550 to 850 
GeV. The range of mass that they found depended on the relative importance of coannihilation 
with the level-one KK excitations of the right-handed leptons. This was the only coannihilation 
channel that was considered in || . Several groups have also examined the prospects for detecting 
KK dark matter directly and indirectly @. 

Coannihilation is likely to occur through other level one KK excitations. Examining a typical 
spectrum of KK excitations from the one- loop radiative corrections discussed in ||, one finds 
that the KK excitations of the particles transforming under the electroweak (but not color) 
groups tend to be within about 10% of the mass of the KK photon. This is near enough in 
mass that coannihilation with these KK excitations is expected to affect the the thermal relic 
abundance. We emphasize that these radiative corrections, proportional to the log of the ratio 
of the cutoff scale A to the compactification scale -R -1 , are merely indicative of shifts in the level 
one KK particle masses, since the finite matching corrections are unknown. Indeed, if the cutoff 
scale were not too large, even the KK excitations that transform under SU(3) C could also play 
an important role in the calculation of the thermal relic abundance. 

What is the cutoff scale of UED models? Generically, the cutoff scale is where the extra 
dimensional theory gets strong. In the four dimensional effective theory this can be estimated 
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by including KK excitations into loops, and one obtains typically AR ~ 10-100. But there have 
also been recent explicit scattering amplitude calculations that suggest the cutoff scale may be 
much lower than previously assumed [jyj]. In any case, it is clear that including the effects of 
coannihilation for the entire spectrum will allow us to calculate the broadest range of KK photon 
masses that lead to a thermal relic abundance consistent with the cosmological measurements. 
We will illustrate the relative importance of different channels by showing the results for specific 
choices of the level one KK spectrum. 

The lower bound on the mass of the KK photon arises by calculating the effects of KK 
particles in loop corrections to precision electroweak data. This was first done by || who found 
a lower bound of about 300 GeV. As we will see, this bound is not saturated no matter what 
spectrum of level one KK excitations one takes, and so we do not need to consider it further. 

In general, the accuracy we desire requires only tree-level calculations. However a very 
recent paper did find an important effect on the relic abundance of the KK photon The 
loop suppressed operators on the orbifold fixed points lead to couplings between pairs of level- 
one excitations and SM particles through resonances of level two (and higher) KK excitations. 
The resonance behavior can roughly compensate for the loop suppression, resulting in a large 



contribution to the total cross section. Ref. [12] investigated the effect of the annihilation 
syi 1 )^ 1 ) _> it through an s channel h^ 2 \ and found that including this process leads to a 
significant decrease in the relic abundance of . Including these loop-suppressed effects would 
make our calculations prohibitively complicated, so in this paper we neglect them. Instead, our 
motivation is orthogonal, namely to understand the implications of including coannihilation of 
all level one KK excitations on the relic abundance. Ultimately, to obtain the most accurate 
estimate of the relic abundance of the KK photon would require taking into account both the 
resonance effect as well as the coannihilation effects that we report in this paper. 



3 The Relic Density of the KK photon 
3.1 Annihilation and freeze-out 

We now briefly review the standard method for calculating the cosmological relic density of a 
stable particle [O, 14]. The number density of a particle in the expanding universe is described 



by the Boltzmann equation 



^"^ + SHn^ = - {a A v) (n^) 2 - (n!f) 2 



dt 

where {o~av) is the thermal average of the total cross section for annihilation of the particle ip 
(assumed here to be equal to the total cross section for producing -0) times the relative velocity 
of two particles, and H is the Hubble expansion rate of the universe. 

We solve the Boltzmann equation to obtain the number density of a massive particle at late 
times. At early times (T S> m^, n e ^ ~ T 3 ), n e ^ is large, and any deviation of n from n eq rapidly 
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goes to 0. At late times (T <C m^, ~ g(m^T/2ir) 3 / 2 e~ m ^/ T ), n eq is tiny since the particles 
are now too massive to be thermally produced. As decreases due to the expansion of the 
universe, the particles eventually become too dilute to annihilate and freeze out at a constant 
density per comoving volume. This freeze-out occurs roughly when the Hubble expansion rate 
overtakes the rate at which the particles annihilate, V = (a^v) ~ H. 

As we will verify later, KK dark matter freezes out at Tf <C ttikk, so that at freeze out 
/ 2 \ 3 / 2 

n eq = g I J e x , where x = ^ j an d g is the number of degrees of freedom of the annihilating 

particle. In this case, (||) cannot be solved analytically. Instead, the Boltzmann equation must 
be solved either numerically, or by means of a standard approximation which gives solutions 
consistent with numerical results to within 10%. Here we shall use the latter approach. 

To solve the Boltzmann equation approximately, we begin by making a simplifying change 
of coordinates. If we assume the universe expands adiabatically, then sa 3 = constant, where s is 
the comoving entropy density given by s = ^-g*T 3 and g* is the effective number of relativistic 
degrees of freedom. In terms of the variables Y = n^/s, Y eq = n eq /s, and x = m/T, (|l]) can be 
rewritten as 

sy = _m {y , _ {Y „ f] (2) 

ax Hx 

where we assume that freeze-out occurs in a radiation dominated epoch, so that 4? = Hx. 
Writing A = Y — Y eq , this becomes 

A' + Y eql = -C^^-A(2Y eq + A) 

= -f(x)A(2Y eq + A) (3) 

where C = J ^mMpj, and f(x) = C (ctav) x~ 2 . 

At early times, the particles track their equilibrium values, so that A <C Y eq and A' <C Y eq/ . 
In this case 

yeql 

A = ~ f(x)(2Y eq + A) (4) 

To determine the freeze-out temperature, we use this expression for A to solve A = cY eq with 
Y eg = ^9 x 3 / 2 e~ x , where c is a numerical factor whose optimal value is determined by 

v327r'g* 

comparison with numerical solutions of the Boltzmann equation (c = v2 — 1). The freeze-out 
temperature can be deduced by solving the following equation numerically: 



At late times, Y eq and Y eql are both small, and the equation reduces to 

A 7 



^2 = -/(*) ( 6 ) 
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This equation is separable, and can thus be solved (with boundary conditions again chosen to 
give the best fit to numerical simulations): 



Y k A 

- 1 no ~ *-* r 



Pi, s ^/45 



Pc ^/KgZM-pipc J^° f (cr A v) x 2 dx 

1.04 x 10 9 



f™ (<TAV) X' 



2 dx 



(7) 



Note that because of the factor of sq appearing in the numerator, this solution assumes that the 
universe expands adiabatically between freeze-out and today. This is, of course, not quite true, 
as some standard model particles will subsequently fall out of equilibrium, slightly increasing 
the entropy per comoving volume. In practice, however, this contribution from known standard 
model sources is very small, and hence may safely be neglected. 



3.2 Coannihilation 



Different species of interacting particles that have masses nearly degenerate with the KK photon 
will fall out of equilibrium at the same time and with roughly similar density. Interactions of 
the form tpiX ^ ip2X' (where if)% and ip2 axe the KK particles, and X, X' are SM particles), 
converting the relic into other particles of similar mass, occur rapidly. In fact, such processes 
are much more efficient than the relevant annihilation processes, so that the abundances of the 
all particles are effectively correlated during freeze-out. The slightly heavier particle species 
will eventually all decay to the lightest stable particle (in our case, the KK photon) and thus 
contribute to the relic density. This situation is called coannihilation. 



As shown in [14], to compute the relic density and freeze-out temperature with coannihilation, 
it suffices to substitute an effective cross section cr e g for a in Eqs. (0),((5|) above, where 

g cS = Y,9 i a + ^f /2 e- A ' x (8) 



The relic abundance depends on a weighted average of the annihilation and coannihilation cross 
sections of all relevant particles. Thus if the coannihilating particles interact more strongly 
compared with the annihilating particles (while also being close enough in mass) they can 
increase the effective cross section and thereby decrease the relic density. Conversely, if they 
are more weakly interacting, then the coannihilation effects decrease the effective cross section, 
yielding a larger final value for f2^. 
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A simple intuitive picture of the effects of including coannihilation comes from noting that 

P ~ ^jT° M x'A . (9) 

The integral over x will of course result in a different weight for s and p wave annihilation, 
nevertheless it is easiest to simply ignore this difference and approximate ~ ((av))^ 1 . A 
larger effective cross section results in more annihilation during freeze-out, and thus a smaller 
relic density. 

Furthermore, it is easy to see that the effective cross section a e g is a weighted average of 
the relevant annihilation and coannihilation cross sections by considering the case where all but 
the lightest (first) coannihilating particles have the same mass mj = mo(l + 5), j 7^ 1, and all 
particles have the same number of degrees of freedom. Then 



^(1 + <J) 3 (^+ j)/ 2 e - ( »+ j)& B £ii 



£(1 + 5f^l 2 e 



- (10> 

v 

where i,j run over all included particles, a^j = 0(1) for the LKP (all other KK particles), and 
Gij is the total (co) annihilation cross section for the process ipiipj — > Standard Model particles. 

This makes it clear that in the limit of a small mass difference and equal numbers of degrees 
of freedom, cr e ff is precisely the average of all relevant cross sections while the relic abundance 
depends only on whether this average is larger or smaller than an. If all cry are equal, then 
including coannihilation does not alter the relic density; if the weighted cross sections of the 
coannihilating particles are on average larger (smaller) than that of the original particle, then 
the relic density will be smaller (larger) when coannihilation is included. This holds irrespective 
of the value of 8. However, since 5 enters exponentially into the weights of the coannihilation 
cross sections, the effects of coannihilation are rapidly suppressed as the relative mass difference 
5 increases. Typically for weakly interacting particle dark matter, xj ~ 20 to 30, and for 5 > 0.1 
the effects of coannihilation are generally negligible. However, for very large coannihilating cross 



sections, fractional mass differences up to 5 = 0.2 can affect the relic density [14]. 



4 Results 



Now we are ready to consider the effect of coannihilation on the relic abundance of the the 
KK photon. Each coannihilation scenario considered requires several total annihilation cross 
sections; the relevant Feynman diagrams, together with tables of the formulae of the relevant 
cross sections, are presented in Appendices ^| and O 

Certain approximations were made in obtaining our results. Cross sections were computed 
ignoring all terms of order !" SM where msM is any standard model mass. Consequently, we 
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ignored all Yukawa couplings which are proportional to the corresponding fermion mass. This 
approximation is well justified (in that including such terms should alter the relic abundance by 
less than 10%) for all particles except possibly the top quark, whose mass mt ~ 175 GeV is close 
to one half of the precision electroweak lower bound on In particular, as shown in [[T^], 

the top Yukawa coupling can alter the cross section significantly, as it leads to nearly resonant 
s-channel diagrams. Aside from this effect, our results are not expected to be sensitive to the 
top mass since, as we will see, the masses of the KK photon necessary to explain the observed 
relic dark matter abundance are well above the mass of the top quark. It addition, we neglected 
the mixing between and VF 3 ^), which are expected to be rather small already for the first 
KK level (see Appendix [O]). Hence, we take 7 (1) = , = W 3 (1) throughout. 

When examining the relic density as a function of mass, a further simplification is obtained 
by ignoring the mass-dependence of xt. The value of Xf depends very weakly on the mass. 
Indeed, from (||) we expect this dependence to be approximately logarithmic. Typically, over 
the mass range of tukk = 0.2 to 2 TeV, Xf varies by about 2 GeV/degree, or less than 10%.* 
This variation has a small effect on the relic density. This also shows that the KK dark matter 
is cold: for all cases considered here, we find 23 < xj < 28, so that the particles are well 
approximated as non-relativistic. 



We use couplings and astrophysical parameters that can be obtained from 15]. In particular, 
we take the electroweak and strong couplings evaluated at Mz (ignoring the renormalization 
group running of the couplings). In practice, the freeze-out temperature is between about 25 
and 100 GeV, depending on which particles coannihilate. The dark matter relic density is taken 
to be the la value from WMAP, n DM = 0.113 ± 0.01 M. 



4.1 Coannihilation with simplified level-one KK spectrum 

To investigate the effects of KK particles coannihilating with the KK photon, we chose two 
different mass spectra for the coannihilating particles. First, particles in the first KK level are 
divided into two groups: those with masses approximately m (i)(l + 5), with 5 small, and those 
that are too heavy to play a role in coannihilation. While this choice of mass spectrum gives a 
good understanding of the range of possible outcomes of including coannihilating particles, it is 
also somewhat arbitrary. Thus, we also investigate a mass spectrum based on that derived in 
H in the next subsection. 

The first check on our results was to consider coannihilation of the KK photon with just 
right-handed leptons, where our calculations and numerical results agree with §. This provides 
a non-trivial check on our procedure for calculating the relic density in the Universal Extra 
Dimensions model. 

Figure |l| shows the effect of including coannihilation of the KK photon with leptons (blue) , 
and all electroweak particles including leptons, scalars, and electroweak gauge bosons (red). 

*The dependence of Xf on the relative masses of the coannihilating particles, however, has a somewhat larger 
effect, and have been taken into account. 
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Here we assume in each case that particles included in the coannihilation have the same mass 
m 7 (l + 5). The graph shows that coannihilation with leptons, and to a lesser extent with W 
and Z bosons, tends to increase the relic density for a given mass of the KK photon. 

Figure |2| shows the effect of including coannihilation of the KK photon with all level one 
KK electroweak particles and level one KK quarks (blue), and all particles at KK level one 
including the KK gluon (red). Again, we assume in each case that all particles included in the 
coannihilation have the same mass m 7 (l + 5). Clearly, including coannihilation with strongly 
interacting KK particles decreases the relic abundance for a fixed KK photon mass. The extent 
of this decrease is one of the most important results of this paper. In particular, if the level 
one KK particles are highly degenerate to within a few to perhaps 10%, the KK photon mass 
consistent with the thermal relic abundance that matches the WMAP data could be several 
TeV. 

These results can be more readily understood by examining the relative magnitudes of the 
cross sections in question, shown in in Fig. ||. Notice first that the coannihilation cross section 
of any non-colored particle with B^ is smaller than the annihilation cross section of B^ with 
itself. Thus coannihilation with particles whose self-annihilation cross sections are close to those 
of B^ tends to decrease the effective cross section, and hence increase the relic abundance for 
a given B^ mass. This is the case for leptons and scalars. The level one KK electroweak 
gauge bosons and and all strongly interacting particles have sufficiently large self- 

annihilation cross sections that they cause an increase in the effective cross section, and decrease 
the relic abundance. For and Z^ this effect is relatively small, and their inclusion does not 
completely counter-balance the reduction in <r e // resulting from including leptons and scalars. 
In the case of the strongly interacting particles, the relevant cross sections are approximately an 
order of magnitude larger than those of the non-strongly interacting particles, and the resulting 
reduction in the relic density is quite dramatic at small 5. 



4.2 Coannihilation with the Cheng, Matchev, Schmaltz spectrum 

In this section we consider the spectrum that results if one takes the one-loop radiative cor- 
rections to the KK masses from ||. There are several assumptions built into this spectrum. 
One is that the matching contributions to the brane-localized kinetic terms are assumed to be 
zero when evaluated at the cutoff scale. Furthermore, the radiatively generated terms are log 
enhanced by a log of the ratio of the cutoff scale A to the mass of the KK excitation. This 



ratio is also not known, and may be much smaller than previously estimated [11 1. We therefore 
consider a wide range of AR to illustrate the effects of coannihilation. 

First, let us summarize the results of M that are used for the masses of the KK excitations 
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Figure 1: Relic abundance of KK dark matter as a function of mass after including no 
coannihilation (black; right-most solid line), coannihilation of with all leptons (blue; left- 
most pair of solid and dashed lines), and all electroweak particles (red; middle pair of solid 
and dashed lines), assuming in each case that all coannihilating particles have the same mass 
m ~(i)(l + $)■ The solid (dashed) lines show the values for 5 = 0.01 (0.05) for the cases with 
coannihilation. Notice that for the case including KK electroweak gauge bosons, the abundance 
as a function of mass is smaller for the smaller mass splitting S = 0.01, unlike the case with 
just KK leptons, since the effects of coannihilation of these two sets of KK particles somewhat 
compensate for each other in the relic abundance calculation. 



10 




Figure 2: Relic abundance of KK dark matter as a function of mass after including coannihilation 
of flW with all non-strongly interacting particles and quarks (blue), and all level one particles 
(red) assuming in each case that all coannihilating particles have the same mass m (i)(l + 5). 
The solid and dashed lines show the values for 5 = 0.01, and 5 = 0.05. 
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Figure 3: Magnitudes of all cross sections, using txikk = 1 TeV, xj = 25, and all couplings 
taken at the scale Mz- Here each color represents a particular particle species (see legend), and 
the x-axis indicates with which particle the (co) annihilation is occurring. For example, a blue 
dot above the text ex is the coannihilation cross section of = 7W with e^. Only 2 families 
of fermions are shown, as this is sufficient to show the difference between annihilation between 
members of the same family, and coannihilation between members of different families. The 
annihilation cross sections for fermions, W and its scalar counterpart w have been weighted by 
a factor of 1/2 to account for the two helicities. In general the coannihilation cross sections 
are smaller than the annihilation ones, (and more numerous) so that usually adding more 
coannihilating particles decreases the effective cross section. However, as the SU(2) coupling 
is stronger than the U(l) coupling, and as it also opens more channels, when scalars and gauge 
bosons, whose annihilation cross sections are quite large, are added, the average cross section 
increases. 
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All of the non-colored KK excitation masses are within about 10% of m (1) up to moderately high 
values of the cutoff scale (Ai? ~ 30), and hence are almost certainly relevant for coannihilation. 
Conversely, for Ai? > 5, the masses of the strongly interacting particles are more than 10% 
greater than m^i), and thus are less likely to be important for coannihilation. However, as em- 
phasized by |lij| ) particles with sufficiently large cross sections may be relevant to coannihilation 
even for mass differences of up to of order 20%. 



We further simplify the mass spectrum given by (11), while ignoring the SM top quark mass, 
by dividing the KK particles into five mass classes: the KK photon; the KK leptons and scalars; 
the KK W and Z gauge bosons; the KK quarks; and the KK gluon. (In cases with more than 
one particle of slightly differing mass, we take the average mass of all of the particles in that 
class.) 

Figure || shows the mass that results in a thermal relic abundance consistent with 
cosmological data, as a function of Ai?, given a five class mass spectrum. The upper plot also 
shows the mass gap 5 for each of the classes. Because g s is considerably larger than g,g', the 
masses of the strongly interacting particles increase rather quickly as the cutoff increases. For 
Ai? > 2, nearly all strongly interacting particles have 5 > 0.05. Masses of the leptons and 
scalars, however, vary rather slowly with Ai?, and remain within 10% of the mass out to 
values Ai? ~ 10 4 . Thus as the cutoff is increased, the gluons, quarks, and to a lesser extent the 
W and Z bosons, rapidly become too heavy to play a role in coannihilation, and the relic density 
is determined by the effective cross section of 7W and the KK leptons and scalars. Thus the 
mass spectrum, Eq. (JTT|) , favors a somewhat lower value of m (1) as 700 GeV to be consistent 
with the measurements of the dark matter abundance in the Universe. 
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Figure 4: The calculated mass yielding a thermal relic abundance consistent with WMAP 
observations is shown in the bottom graph using as a function of AR. This result was 
calculated using a simplified spectrum of the first KK level consisting of four sub-levels that 
vary logarithmically with AR, shown in the top graph. The ratios are from top to bottom: the 
KK gluon, the KK quarks, the KK electroweak gauge bosons, and the KK leptons and scalars, 
divided by the KK photon mass. The black dashed line corresponds to 5 = 0.05 for reference. 
In the bottom graph, including increasing numbers of particles at the first KK level are shown 
with the four lines labeled in the legend. The solid straight black line (the line at the top on 
the far right) shows the the case without coannihilation for reference. All level one KK particles 
except for the leptons and scalars become rapidly irrelevant to the relic abundance calculation 
once Ai? > 100. 
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5 Conclusions 



We have calculated the thermal relic abundance of the KK photon in the five dimensional UED 
model for a generalized mass spectrum of level one KK particles. We find that the lowest KK 
photon mass which could possibly account for the observed dark matter relic abundance is 540 
GeV, resulting from including coannihilation with three generations 6r and neutrinos all very 
nearly degenerate with 7W (# = 0.01). This is consistent with the result found in Ref. ||] 
updated to reflect the WMAP measurements. Since the radiative mass corrections should be 
the same for ei and v, a more realistic estimate, given by including all KK leptons, is 570 GeV. 
This is significantly lower than the lower mass bound of 860 GeV given by including alone. 
Including level one KK quarks and the KK gluon with masses within about 10% of the mass of 
the KK photon significantly increases the total effective annihilation cross section. This implies 
the KK photon mass leading to a thermal relic abundance consistent with WMAP observations 
is much larger, up to several TeV, see Figs. || and || for more precise numerical ranges. On face 
value, such a small separation between the KK photon and the strongly interacting level one 
KK particles is not expected from the radiative corrections to the masses of the first KK level 
computed in p|. However, if the cutoff scale is not much larger than the KK photon mass itself, 
and thus matching corrections are comparable in size while opposite in sign to compensate, the 
level one KK spectrum could be much more degenerate. These results show that the range of the 
KK photon mass is much wider if indeed the mass spectrum is more degenerate than previously 
thought. Given measurements of the level one KK particle masses, the calculations presented 
here could be used to find the total effective cross section to verify if the KK photon does (or 
does not) make up the dark matter density needed to be consistent with WMAP observations. 

Note added: During the completion of this work we became aware of an analogous 
calculation done independently by another group ||16|| . We have compared extensively the 
formulas for the annihilation cross-sections involving KK quarks, KK leptons and KK gauge 
bosons, in the limit of degenerate KK masses, as listed in the Appendix. In all considered cases 
we found perfect agreement. 
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A Universal Extra Dimensions 

Here we provide a brief overview of the Universal Extra Dimension model in five dimensions (for 
a review, see [|l7j). In subsequent Appendices we list the Feynman diagrams for (co) annihilation 
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and the cross section results we obtained. The particle content of the UED model is shown in 
Table [l]. The Z2 orbifold projects out one helicity of the fermion zero modes, as well as the 
zero mode of the 5 th component of the gauge field. We begin with a brief survey of the mass 
eigenstates of the theory; the Feynman rules are presented in Appendix M. 



5D field 


AD Even Fields 


4D Odd Fields 


A M 




A 5 




L = (u, ei) 


L = (vr, e R ) 


(Q,Q) 


Q = (u L ,d L ) 


Q = {u R ,dR) 


(e,e) 


e = e R 


e = e L 


(u,u) 


U = Ur 


u = u L 


(d,ij 




d = dR 

4> 


d = d L 



Table 1: 5D fields and their behavior under the orbifold projection y — > —y. 



A.l Fermions 



The kinetic terms for fermions have the form: 

Luetic = ip(x, y)(7^ + 55,57555)^(2;, y) (12) 

where we have ignored SM mass terms. Upon doing the usual Fourier expansion for, for example, 
the first generation of (left-handed) quarks, and integrating over y, we obtain: 



Luetic = (u, d)L7 M 5 M (n, d) L + £ P L Qf L h^P L Q[ j l - J^PrQH 

j= l L V 



+ p RQm {fd*PRQxR + ~1% p LQil) \ ■ ( 13 ) 

(We use lower case letters to denote SM particles, and upper case for their KK excitations.) 
To obtain the full fermion masses, the usual mass terms arising from the Yukawa couplings are 
added (then diagonalizing the resulting mass matrices). 



A. 2 Gauge Bosons 

After compactification, the kinetic terms in the gauge boson Lagrangian can be expressed as 

L kin etic = -~ {F^F^ + 2{d b A lx - d^a) 2 ) (14) 

where A^ is the 4D gauge field, and (to make the analogy with the Higgs mechanism more 
apparent) a = A§, a 4D scalar. Integrating over the fifth dimension coordinate y causes all 
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cross-terms between modes of different KK-number to cancel, leaving 

oo oo ^ . ^2 

W-tJ 8^ 1 -H W +oE • ( 15 ) 

3=0 3=1 V J 

Expanding the second term, and collecting modes of different KK number, we obtain: 

3=1 



1 2 1 

+ - ( 41 A&'U^ - f-1 9 tf a^U"^ + -duaWdfa® 
2\RJ 11 \RJ 2 



(16) 



This is precisely the Lagrangian for a tower of 4D gauge fields with masses Mj = ^ generated 

by a spontaneously broken symmetry. The scalar fields a® are eaten by the gauge field in 
the usual Higgs mechanism to give them mass. 

A. 3 Electroweak KK Gauge Boson Mass Eigenstates 

There is a tower of KK gauge bosons for each of the SM gauge symmetries. The electroweak 
gauge boson KK tower, however, is different from the SM zero modes in the admixture of W 3 C?) 
and BV\ This is because the KK SU(2) and U(l) gauge bosons receive different loop corrections 
to their masses. The actual KK mass eigenstates can be found by diagonalizing the matrix 



+ 6i + \g 2 v 2 \gg'v 2 



(17) 



where 8\ = 8{m w ( n )) 1 82 = 8(m B ( n )) are the radiative corrections to the n th KK level electroweak 
gauge bosons (precise expressions can be found in [f|). We can re-express this as: 



,2 



v £ + a ai j T V ba 1 " 1 s + is"'"- 

where Am = Sl ~^ S2 , and 8 = Sl ~ s * . The first contribution to the masses is diagonal in the 
W^-B basis, and does not affect mixing. The mass eigenstates result from diagonalizing the 
second matrix; thus the Weinberg angle at each KK level is determined by the relative sizes of 
8 and the electroweak masses \gigjV 2 . Since 8 is proportional to if thkk ^> M ew , the KK 
mass matrix is almost diagonal in the W^-B basis. Since tukk is at least 300 GeV Q, the KK 
Weinberg angle is much smaller than its SM counterpart. In fact, for R^ 1 > 600 GeV (which 
is roughly the lower bound on the mass of B^) and A/fi > 2, sm(6 KK ) < 0.15, and to the 
accuracy required here we can take ~ B^\ and ~ W^ 1 ^- 
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A. 4 Scalars in the Electroweak Sector 



a (n) 

We have just shown that in the absence of other scalar interactions A 5 is the Goldstone boson 
eaten by in the effective 4D theory. In the weak sector, however, the interactions between 
A$ and the KK scalar fields mix up the mass eigenstates. Note that this was also discussed in 
Ref. PL 



To see this, consider the 2-point vertices of the SU(2) gauge bosons. The 5D gauge field 
strength contributes a term of the form (|l6|). Now we add to this the interactions with the 
standard model scalar fields: 



L 



scalar 



\\dm<S> 



igA a m T a <5>\ z 

±{d m <S>d m $ + igA a m [(T a ^d m $ - d m &T a $ 
+ g 2 A a m {T a ^A^{T b $>)} 



(19) 



where A m , $ are the five-dimensional fields. We now will determine the physical and Goldstone 
scalars by examining all of the 2-point vertices in this Lagrangian. 

After integrating out the 5 th dimension, and adding in the gauge boson kinetic term, the 
relevant piece of the Lagrangian is: 

mM4a^4 n) - -\<S>\ 2 m KK 2 + -A$ rr<W$ (n) * - d^^ n) {T a v) 



+ - 



+| a <W \ T a vm KK & n > -m KK & n \T a v)* 
AtfriA™ \(T a v)^T b v + m 2 KK \ - a a (n)a^ (T a v^T b v 



(20) 



where v = (0, v) is the standard model Higgs VEV. The terms involving T a vdM& give 2- 
point couplings between the KK weak gauge boson and the KK excitation of the corresponding 
Goldstone boson, which we write as 3>i"^. In the case of A|(= a a ), these 2-point couplings 
involve p§ = ttikk rather than spatial derivatives, and thus become mass mixing terms between 
$ a and a a . In addition, both $ and A® acquire KK mass from their interactions with a a , while 
a a gets a mass from its interactions with the Higgs VEV. Thus a a is no longer purely a Goldstone 
boson, as it has an electroweak scale mass. 



A convenient gauge- fixing functional is: 



G 



1 



d^A^-i{gF a ^+m KK a^ 



(21) 



After fixing the gauge via L — > L — \G 2 , all 2-point interactions between the scalars and gauge 
bosons of the effective 4D theory are canceled for any gauge fixing parameter £. This means 
that one is free to choose a gauge (£ — > 00) in which the Goldstone boson propagator vanishes. 
For general values of £, there are mass mixing terms between a a n ^ and This leads to mass 
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matrices for the scalars of the form 

m2 KK^ + m z -mzmKK (£ - 1) 
-m z m KK (£ - 1) m|£ + m\ K 



(22) 



between a a and the KK excitation of the corresponding SM Goldstone boson. Here m\y and 
mz are the masses of the SM W and Z gauge bosons, respectively. 

To obtain the physical and Goldstone particles, we diagonalize the mass matrix.''' We find 
the mass eigenstates 

p^ = (—rriKKd + m z (fr) with mass £ yj m? KK + m 2 z 

\J m2 z+ m KK 

p^ = (m z a + rriKK^) with mass \J rn? KK + m? z 

I ml + m 2 KK 

Clearly pW i s 

the Goldstone boson, which we can be eliminated from the theory by taking the 
limit £ — > oo, leaving p^ the physical scalar particle. The fact that the physical scalar is a 
mixture of A^ a and $ a changes its couplings to other particles by terms proportional to the 
masses of the electroweak gauge bosons. Note that this mixing in the scalars is needed to ensure 
unitarity of gauge boson scattering is preserved, which we explicitly verified. 

Since the A^ component of the physical scalar is suppressed by relative to the 

component, in practice this mixing can be ignored in all calculations not involving massive 
external SM gauge bosons. 

This completes our discussion of the particle spectrum of the UED model. The KK zero 
modes are the particles of the Standard Model. At ukk = 1 and higher, the effective theory 
contains massive vector bosons that have eaten the corresponding Goldstone (or a combi- 
nation of A^ and for the W and Z). It also contains both helicities of SU(2) doublet and 
singlet fermions, with the SM helicity even under the Z2 action, and the non-SM helicity odd. 
Finally, the model contains physical scalars, which are the KK Higgs, M n \ and that 
are mixtures of A^ ,Z with the KK excitations of the SM Goldstone bosons (f>w,z- 



B Feynman Rules 

In this section the Feynman rules relevant for the calculations used in this paper are written. 
Only the relevant interactions, namely between SM particles and level one KK modes, are 

^To avoid potential confusion with signs, note that the signs of the off-diagonal terms in the mass matrix depend 
on how we define $ and For the above, we have used the parameterization $ = (i(wi + iw2),v + <f> — iz); 
interchanging "!> and $* in this convention will change the sign of the off-diagonal terms in the mass matrix, which 
introduces a relative — sign between the a and $ components of the physical scalar. In order to ensure that the 
cancellations necessary to maintain the unitarity of the theory occur, care must be taken to use the correct sign 
conventions when computing vertex couplings. 
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shown. In the diagrams we use double and single lines to denote KK particles and SM particles, 
respectively. 



B.l Fermion/Gauge Boson Interactions 



The fermion interactions of the KK modes differ from those of the standard model due to the 
vector-like nature of these higher modes, which introduces helicity operators at certain vertices. 
We find it convenient to work in unitary gauge, so that the Goldstone bosons do not appear as 
external particles but instead as the longitudinal polarizations of the massive KK gauge bosons. 
The relevant fermion interactions, after integrating over the fifth dimension, are 



/ttR 



fW-fAWp L Fjp + P^if V4 X) / (0) 



W 

/ttR 



P L F^AfP L F^ + P R F^AfP R F$ 



leading to the following vertices: 
f(o) F (D A W 




(23) 
(24) 



igj^Pi (ingoing / or outgoing F) 
igj^Pji (ingoing F or outgoing /) 



F^F^Af 




where FlPlFl + FrPrFr = FF; i.e., this vertex is non-chiral as expected. 



B.l.l Fermion/ Scalar Interactions 



The interactions of fermions with scalars proportional to Yukawa couplings are presented for 
completeness, even though we do not make use of them since we have ignored terms of order 
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v 2 /m 2 KK . The Yukawa terms result in 



7 (0) 



K ,P R G^ + ^=P R G^ 



P R F R l) 



'■kR 
K 
/ttR 



I-kR 



PlF? 



/ttR 



$*(!) 



+ 



nR 
Xd 
V^R 



9d 



(25) 
(26) 
(27) 
(28) 



where F (/) and G (g) are the KK (SM) SU(2) singlet and doublet fermions, respectively. In 
addition, the interactions of A§ with fermions can be deduced from the gauge boson/fermion 
interaction Lagrangian above. The only difference is that A§ is odd under the orbifold Z 2 , and 
so it couples to the "wrong" (opposite) handedness of the KK fermions. After integrating over 
the fifth dimension, the interaction terms become: 



/ttR 



and similarly for G, replacing left with right KK fermions. This leads to the Feynman rules: 
f(o) F (i) A W 

^P R (for F) 
Pl (for G) 



(29) 




/tt R 



/(o)<2(i)$(i) 




/ttR 
Ar 



i?(i) a (o)$(i) 



^^Pr (for G u ,$*) 
^(for G D ,$) 

ia 2 P L (for g u ,$*) 
B=P L (for g D ,$) 



'tt R 



/tt R 



^(1)^(1)0(0) 



^ia 2 (for G u: . 
(for G D , 4 



1 7T R 
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where FlPrGr + FrPlGl = FG. Unlike the vertices between SM scalars and fermions, which 
are suppressed by the small Yukawa couplings, the intrinsic strength of the vertex between A5 
and fermions is not small. However, since the physical scalar particle is mzAr, + itikk®, in the 
limit that that both fermion masses and mz are small compared with ttikk, we can ignore all 
scalar- fermion vertices. 



B.1.2 Gauge Boson/Scalar Interactions 



The electroweak gauge boson/scalar interactions are of the form 



+ 



9_ 

ttR L 



AWA^ + $*0 + v<5> + v<5>*) + ij, 1 ^' 1 '" (\<P\ 2 + W> + vcjf) + ^°)^°^|$| 2 1 (30) 



where the gauge group indices and generators have been suppressed. The form of these interac- 
tions are nevertheless identical to those in the Standard Model. 



Since A§ and the Higgs mix with each other, as we showed above, it is convenient to compute 
the interactions of ^5 with electroweak gauge bosons and scalars. In the electroweak sector, 
where f abc = ie abc , these interactions arise from 

ig(ff>Ai - d 5 A^)(AlAi - A b 5 Al) - ^{A\A% - A\A^ . (31) 

After integrating over the fifth dimension, the only terms which survive are (using the notation 
afc to denote A^ k ): 

(-d»Af n) m { $ K A a ^ (A b ^Af n) - Af l) Af">) - ^ {A b ^ A% n) - A b ^ n) A c ^f (32) 

Summing over a,b,c € SU(2) and switching into the usual W + ,W~ , A, Z basis gives the 
following 3-point interaction terms: 

i (- W +i - n \W- v + W-( n K W + v ) m KK (s w a 7 + c w a z ) 
+i(w+W v {Z vCw + A v s w ) - {A { l l) s w + Zl n) c w ) W+ v ) m KK a_ 

+i (- W^ n K(Z u c w + A v s w ) + (Zl) n) c w + ^ n) s w ) W-J) m KK a+ (33) 

for the two gauge boson/^ interactions, and 

-i (a^daf - a?W n) ) s w A^ - i (-af } £>a_ + a^da^) c w Z^ 
-i ({ai n) c w + a^s w ] da^ ] - aL n) [da^ s w daf> 'c,}) W+(°K - 
1 (- {ai n) c w + a^s w ) da™ + a? [da™ s w + da^c w }) W~(°K (34) 
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for the triple-scalar /gauge boson interactions. The 4-point interactions are obtained from the 
corresponding standard model vertices by requiring both KK particles be A§. 

A§ also couples to $ through |Z^ M <I>| 2 . This results in three types of vertices: V\ = grriKK A5 
V2 = gmzA^A^, and V3 = g 2 A5A5 <£<!>*. V2 and V3 are identical to those of the KK gauge 
bosons, once the factors of g^ u are replaced with #55. V\ is obtained from the corresponding KK 
gauge boson vertex by replacing — > itxikk- As some care must be taken with signs, the V\ 
terms are listed here: 

- is w (-w+Ww- + wJ n ^w + ^ a$f> 

+ ^7r^4 n) ( i cos (26 w ) (w+ {n) u»_ -wJ n) w + ) + h^z-z^h) 

2 cos (p w ) V / 

+i (izu;J n) - hwJ n) + /i (n) w_ - ^u>_) a + (35) 



The vertices involving physical scalars follow from combining vertices involving A§ and in 
the correct proportion. A table of these physical vertices is given below. In principle, care must 
be taken in vertices involving the Weinberg angle, as the mass mixing is not the same for the 
KK particles as it is for the SM particles. In the case of the scalar A$, however, the physical 
component of A$ is that which mixes with the SM Goldstone boson of the Z particle, and hence 
its mixing angle should be the same as that of the SM Z boson. For the KK gauge bosons, the 
mixing angle is different. However, as we already discussed above, it is a good approximation 
to neglect this mixing and thus take ~ and Z^ ~ VF 3 W. 



The following notation is used in the table. We define k w = , KK and k z = , KK . 

^mj iK +mf v V m KK+ m z 

$2,$+, and <&_ are the physical level one KK scalars that are themselves mixtures of the KK 
excitations of the SM Goldstone bosons and A$ from the weak gauge bosons, as described in 



Sec. A. 4 . (Here we use level one KK particles, but the vertices are the same for level n KK 
states also.). §h is the level one KK Higgs particle, while h is the SM Higgs particle. The table 
includes only vertices which differ from those of the SM. 
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2 vector /scalar vertices 




Coupling 

-iCiui-lfrrizKzg^ 
ic w (-lfm z K z g^ u 
im z (s 2 w + {-If 
-im z (s 2 u + (-lfc^Kwg^ 
-im z {s 2 w - {-i fcD^g^ 
im z {s 2 w - {-\Yc 2 w )K w g^ v 
\im z sva.(26 w )((-lf - l)K w g^ v 
^-imzsin(26» u ,)((-l) /3 - l)n w g^ 
\im z sm(26 w )((-lf + l)n w g^ 
=±im z sm(26 w )((-l) 13 - l)K w g^ u 
as in standard model 



2 scalar /vector vertices 



Vertex 
4 0) * + <l>- 

W^ 0) <t> z <S>, 



g{c 

9 

-9 



Coupling 
gs w (p+ -p-) 

cos (28 w ) 2 \( 
2c w K w)\P-\ 



P- 



+ o K z K, w {p- -p z ) 

Kj x SM vertex 



2 vector /2 scalar vertices 

IT 



< (0) W- (0) $ + $_ 

A^A^°h + <t>_ 

z< 0) w*°><m> ± 

4 0) ^ ±(0) $^ ± 



— U + ^T—) K w9tiu 



314(1 + i^) 9i 



2^2 ^cos(26>^,) 2 



3 «: 



"2^ 



+ 



iiij 



2 C 4 m 2 



2 /£os(2M 
9 eK w\ c„ 



+ 



S 2 



'"KK 



1 — 

~)9, 



)9 



2" )Huv 



3 scalar vertices 



$ z <S> z h 


-(m|(l + 




+ 


W 2 ) 




m KK 
m w 


c w rn z 


<S> w <f> w h 


-{m 2 w {l + 


) + 




K 2 


m KK 





4 scalar vertices 



® z $ z hh 



■ 



P^i^P- + 2^- 

1 V m 2 m 2 k . 

2,„2 I\.f2 2 
( h _l_ O Z 



7/: 



Here we have included only vertices in which the SM scalars are the physical Higgs particle. 
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In practice, to lowest order in v/ttikk, the only vertices that are affected are those whose 
analogue in the SM contains a factor of mz- In this case, the contributions of certain A§ terms 
can be of the same order in mz as the contributions of $>W. 

Note that vertices of the form A^A a c have a sign which depends on whether the indices 
o, b, c occur in clockwise or counter-clockwise order, caused by e abc from the commutators of 
SU{2) generators. Diagrams of the form A a ^A^, have no such sign. Hence care must be taken 
in combining the contributions of these two types of diagrams to obtain the physical scalar 
vertex. In order to make this clear, we have explicitly kept a factor of (— 1)^ in the relevant 
vertices. Here (3 is even for vertices in which the electric charge entering the vertex increases in 
the counter-clockwise direction, and is odd otherwise. 

The gauge boson self-interaction terms are exactly as in the standard model, and need not 
be reviewed here. Note that A^ , being even under orbifold parity, has no vertices with the 
Goldstone boson A§ apart from those in the kinetic terms described above. 

C Generic Diagrams and their annihilation cross sections 

In this section, we list all Feynman diagrams involved in the annihilation processes we consider. 
(As explained above, processes involving the Yukawa coupling of Higgs to fermions, the self- 
coupling of Higgs to itself, are not calculated as they are expected to have a small effect on these 
results.) The possible annihilation processes are classified according to the nature of the initial 
and final state particles. Some diagrams apply to multiple processes, and not all diagrams in a 
given section necessarily apply to all processes of that type. A comprehensive list of processes, 
the relevant diagrams, and the corresponding cross sections can be found in Appendix [D]. 

C.l Fermion Annihilation 




C.1.2 - A^A^ 

Here the set of diagrams depends on the particular final state; i.e., the last diagram is absent 
when the final state gauge bosons are hypercharge. 
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C.2 Gauge Boson Annihilation 
C.2.1 A^aP _> /(0)/(°) 

Here vertices involving both SM and KK fermions are chiral, as previously mentioned, and 
factors of the appropriate helicity projection operators must be included. In the limit that 
electroweak breaking masses are ignored, all of these scattering processes result in two final 
state fermions of the same chirality. 




C.2.2 W^WP -> 0«V (O) and Z<P -> 0«V (O) 

This set of diagrams corresponds to the annihilation of W or Z bosons into at least one physical 
SM Higgs boson. In the limit that the masses of the SM W and Z are small compared to 
s, we can consistently ignore diagrams with vertices of the form WWcj), ZZ<p, and WZ<f>, as 
their couplings are suppressed by v/ttikk or v 2 /s. If one of the external particles is a gauge 
boson, the Goldstone boson equivalence theorem tells us that only the longitudinal polarization 
of the external gauge boson contributes. We can then use the Goldstone boson approximation 
to calculate the relevant cross sections from the set of diagrams below. 
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This set of diagrams corresponds to KK gauge boson annihilation into SM gauge bosons. Note 
that we treated SM electroweak gauge bosons as massive, so that no external Goldstone bosons 
are necessary. However, all diagrams with physical scalar propagators must be included. In 
the s-channel, only the Higgs plays a role. However, since there are physical KK w and z 
scalars, these propagators must be included in the t and u channels. We remark that would-be 
high-energy divergences of the form s/m 2 KK are absent precisely because the the non-Goldstone 

scalars are a mixture of and 




C.3 Higgs Annihilation 
C.3.1 -+ 0(0)0(0) 

This set of diagrams correspond to scalar annihilation into at least one final state Higgs particle. 
If final state gauge bosons are present, the resulting diagrams contain at least one coupling 
proportional to mz, and hence are highly suppressed unless the gauge boson is longitudinally 
polarized. Thus all such processes can be described by purely scalar final states in the limit that 
we ignore terms proportional to mz/niKK- 

Here we furthermore ignore diagrams involving 3-scalar vertices, since their contribution is 
suppressed by factors of v/itikk or v 2 / s. 
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C.4 Fermion/gauge boson coannihilation 



C.4.1 f^AP -> /(°)^ 0) 




C.5 Gauge boson/Higgs coannihilation 
C.5.1 cDd)^ 1 ) _> $(o)4°) 

Here again we work in the limit my/ <c s, tukk- In this limit, processes to two final state scalars 
can be neglected, as they necessarily involve at least one vertex suppressed by a factor of mz- 
In processes to one final state scalar and one final state gauge boson, diagrams in which scalar 
propagators are replaced by gauge boson propagators are also electroweak suppressed (by m 2 z , 
while the longitudinal gauge boson contributes \jraz) and thus neglected. 




C.6 Fermion/Higgs Coannihilation 
C.6.1 $(°)/(°) 

This coannihilation diagram includes processes to external Goldstone bosons which are of course 
just the leading order contribution to the analogous process of coannihilating into external SM 
gauge bosons. 




D Cross Sections 
D.l Fermion Annihilation 

The following sections refer to the diagrams and matrix elements computed in Sec. [C], by the 
Section numbers assigned to each process. The total cross section is the sum of the appropriate 
terms and cross-terms, multiplied by the necessary spin averaging and identical particle factors. 
We use the notation L = In ( l-^jt ) . 
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D.l.l e R 



Process 


Diagrams 


Total cross section 


ErEr — > e#e# 
ErEr -> ee 

ErEr — > // 
->■ e RWh~ 
ErUr — > e#/XR 


C.l.l a,b 
C.l.l c,d 

C.l.l d 
C.l.l c 
C.l.l a 


g 4 y 4 

32^5^— [/?s(2s ™Lr) + ™Lr( 4 s $rn 2 KK )L] 
T g^[/3(12 S 2 + 16m^ + 65 S m^) - 12m^(4 S + 3m^)L] 

„4y2 y2 

1 247r/3s 2 V s + KK) 
9iYi R {s + 2ml K ) 2 2 

^^-[/3(4 S + 9m^)-8m^L] 
6iSlfc( 4 «-3<i,) 


-► ^* 


C.1.3 


2 !2 


^R-Ei? -»■ 77 

ErEr — ► 

-► WW 


C.1.2 a,b 

C.1.2 a,b 
C.1.2 c 


{(s 2 + Am 2 KK s - 8m% K )L - s(3(s + Am 2 KK )} 

4 /2(-l/2+sin(e w ))\ 4 

1 8^ ^ {(^ 2 + ^kkS - 8m% K )L - s(3(s + 4m 2 KK )} 

e 4 (s+2m KK 2 ) 
192tt s 2 f3c w 4 



Note that in the s-channel diagrams the mass of the propagator was ignored (i.e., only 
hypercharge gauge photon exchanged), since we neglected terms of order v 2 /m 2 KK . 



D.1.2 e L 



Process 


Diagrams 


Total cross section 


E L El^e' L T L 


C.l.l c,d 


<?s9t 2 e [5s/3-2(2s+3m KK 2 )L] g 4 e [P (4s+9m KK 2 )-8m KK 2 L} 
327rs 2 /3 2 1 647Tsf3 2 rn 2 KK 
g 4 [s+2m KK 2 } 
1 24vrs 2 /3 


e l ~eZ^w + w~ 


C.1.2 a,c 


4 (96 c w 4 m KK 2 +24 sc w 4 )L+2 f3 m KK 2 +sl3-A0 sc w 4 f3-176 m KK 2 c w 4 


9 768tts 2 P 2 c w 4 



Cross sections identical in nature to those of Er have not been reproduced in this table. 
They can be obtained from the previous table with appropriate replacements of the coupling 
constants. 

In the cross section to leptons of the same generation, explicit coupling factors have been 
kept since the diagrams are the same for both (i) ElEl — ► eifiL and (ii) ElEl — > Ve&ei 
as well as the analogous processes with neutrinos in the initial state (iii) and (iv). Here 
g 2 ^ is the total coupling in the s-(i-)channel. Ignoring the SM gauge boson masses, g 2 = 
g 2 = g 2 (tan 2 9wYe l Y £l + T^{Ei)T^{ei)) for annihilation into the same particle, and g 2 = 
g 2 (tan 2 9wYe l Y Ul + T^{El)T^{i/l)) g 2 = g 2 /2 for processes to the electroweak partner (where 
the ^-channel is mediated by a W ± boson). 
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For cross sections to gauge bosons, the leading terms in the limit of vanishing gauge boson 
mass have been kept. 

D.1.3 quarks 

The processes QQ — ► qq and QQ — ► qq are essentially the same as the corresponding processes for 
electrons shown above, except that an additional diagram involving gluon exchange is present. 
In the limit that all KK gauge bosons have mass, and we ignore the masses of all SM gauge 
bosons in the s channel, the relevant amplitude comes from the electron amplitude (with the 
appropriate modification of the hypercharge value) plus the same amplitude with all couplings 
changed to g s . Annihilation to scalars and electroweak gauge bosons are identical to the electron 
case, with the appropriate couplings. Thus the only novel diagrams are annihilation to gg and 

91- 



Process 


Diagrams 


Total cross section 


QQ 91 
QQ^gZ 

QQ gg 


C.1.2 a,b 
C.1.2 a 
C.1.2 a,b,c 


e 2 g 2 ( y -$,m KK A +Am KK 2 s+s 2 )L~s 2 f3-AI3m KK 2 s 

2 I8ws :i l3 2 
9s9 2 z, q (s 2 +4sm KK 2 -8m KK 4 )L-s 2 l3-4:Sl3m KK 2 

36 I S i TT 

g i -4(m KK ' t +4sm KK 2 +s 2 )L+7f3s 2 +31sf3m KK 2 


216 7rs a /3 2 



Here gz, q is the coupling of the quark to the Z boson. Note that in annihilation to gluons, 
the appropriate color factors have been included, along with a factor of § for the average over 
quark colors. In qq — > 57 or qq — ► gZ, a factor of 1/2 must be included for the trace over t a t b . 

D.2 Gauge Boson Annihilation 

Because the KK Weinberg angle is sufficiently small that its effect on our computations can be 
neglected, the results here are presented in the W 3 ^ basis. 

D.2.1 Hypercharge 



Process 


Diagrams 


Total cross section 


B„B V - // 
B^B,, -► W+W- 
B^By — > ZZ 
B^B V hh 


C.2.1 a,b 
C.2.3 d,e,f 
C.2.3 d,e,f 
C.2.2 a,b,c 


g*Yf[(10m KK '+5s)L-7sl3] 

72s 2 tt(3 2 
1 s w 4 g 4 
192 /3c w 4 tts 

1 1 s w 4 g 4 

2 192 /3c w 4 7ts 

1 1 s w 4 g 4 

2 192 0c w 4 tts 
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D.3 W 3 



Process 


Diagrams 


Total cross section 


W 3 ^W 3u - ff 
W 3 ^W 3u - 00 

w 3M w 3l/ - 


C.2.1 a,b 
C.2.2 a,b,c 
C.2.3 d,e,f 

C.2.3 a,b,c,d,e,f 


g*T 3 (f)*[{10m KK *+5s)L-7sl3\ 
72s 2 nf3 2 

1 9 4 

2 192vrs/3 

1 <? 4 

2 192/3 7ts 

ff 4 (-12 m XK 4 s+24 m*^ 6 )L+4 /3 s 3 +3 /3 m KK 2 s 2 +12 sm KK 4 (3 


18 7r s 3 fPmwK 2 
^192 /3tts 



D.4 W ± 



Process 


Diagrams 


Total cross section 


W+W- - f R f R 

w+w- - 
w+w- - 


C.2.1 c 
C.2.1 a,c 
C.2.3 a,c,d,e,s 

C.2.3 a,b,c 
C.2.3 a,b,c,d,e,f 

C.2.3 a,b,c 
2.2.3 a,b,c 
2.2.3 a,b,d 


to leading order in the SM z boson mass 

g 4 (5s+12m KK 2 )L-10l3s-8f3m KK 2 
8 72vr s 2 f3 2 
g i (-32 m KK 2 s-48m KK 4 )L+16 s 2 (3+41m KK 2 (3 s+87 m KK 4 (3 
144 it s 2 /3 2 m KK 2 

e 4 (-12 m KK 4 s+24 m KK 6 )L+4 /3s 3 +3/3 m KK 2 s 2 +12 m KK 4 sf3 
2 18-k l3 2 m KK 2 s 3 
g 4 c w 4 (24 LniKK^+U m KK 4 s(3-l2 Lm KK 4 s+3 s 2 f3 m KK 2 +4 s 3 /3) 
36-7T s A 2 m,KK 2 

^384 13 its 

e 2 c 2 w g 2 (-12mKK 4 s+24m KK 6 )L+4l3s 3 +3l3rn KK 2 s 2 +12sl3 mKK 4 


18 tt s 3 l3 2 mi(K 2 

1 / 

2 192?rs/3 

576 tts 2 /3 



D.5 Gluon 



Process 


Diagrams 


Total cross section 


// 


C.2.1 a,b,c 
C.2.3 a,b,c,s 


^4 (40s+98m Kif 2 )L-83/3s-72/3m KK 2 
64 6vrs 2 /3 2 
4 (-8m KK 2 s 2 -24m K K 4 s+24:mKK e )L+8s 3 f3+13mKK 2 l3s 2 +34mKK 4 sl3 


9s 32tt s 3 l3 2 m KK 2 



Here the factor of 1/64 was included for the average over initial gluon colors. 
D.6 Scalar Annihilation 

In this set of cross sections, note that the KK modes of the SM Goldstone bosons are in fact 
separate propagating degrees of freedom, and thus their annihilation cross sections must be 
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calculated independently of those of the KK W and Z bosons. For processes to one Higgs and 
one SM gauge bosons, we use the Goldstone boson approximation, as processes to final state 
transverse polarizations are suppressed by v/uikk- 



Scalar Annihilation into Scalars: 



Process 



Diagrams 



Total cross section 



H 2 (4f3m KK 2 -8m KK 2 L)H 8/3s+40m KK 2 

512?r sf3c w 4 " 512tt P 2 sm KK 2 c w 2 h 512?r /3 2 sm KK 2 
1 24 s 2 /3-92 fS srriKK 2 -}^ |3m g 4 
768 7r j3 2 s 2 miiK 2 

-2Am KK 2 sL+Aip smK K 2 +16 m KK 4 , 1_ -4 m KI j 4 +f3 sm KK 2 

7r fi 2 s 2 muK 2 c w 2 768 7r /3 2 s 2 rnj<K 2 c w 4 

1 8s/3+4>3m M 2 -4 j3 m KK 2 H 2 +8m KK 2 H 2 L+H 4 (3m KK 2 



W + W — > /l/l 
— ► z/i 

ZZ — > /l/l 
/l/l — > hh 

hh — > z/t 
/iz — ► /l/l 
w; + z — > 



/iz — > zh 



C.3.1 a,b,d 
C.3.1 a,b,c 

C.3.1 a,b,d 
2.2.3 d 



C.3.1 a,b,c 



C.3.1 b,c,d 



C.3.1 a,c,d 



^768 



512 



512 7rs/3 Cu , 4 
~ / J2lz_\ 



r. 



! g 4 (76sm KK 2 /3-4/3m KK i +60s 2 /3-14Asm KK 2 L) 

768 7r /3 2 s 2 mpcK 2 

x 3 4 (84sm KK 2 L-48s 2 /3-30sm xif 2 /3) 



+ 



768 



it f3 2 s' 2 m K K 2 Cw 2 



j_ g 4 (3 sm^ 2 /3+12 s 2 /3-12 sm M 2 L) 

"'"768 -n f3 2 s 2 m KK 2 c w 4 

- 4 12s 2 / 3-20 / 3sm gK 2 -4/3m gK 4 



H 4 



X 

768 it 2 s 2 nriKK 2 

g 4 (l2 sm KK 2 L-6 /3 sm KK 2 ) H 2 ,j _ _ 
+ 768 nP 2 s 2 m KK 2 c w 2 r 256 7rs/3c m 4 

g 4 (-24 sm K y 2 +12 sm KK 2 H 2 )L 
768 7r(3 2 s 2 m KK 2 c„ 4 



+ 



g 4 3 ff 4 /3 sm KK 2 -6 sm KK 2 l3H 2 -4l3m K K 4 +28 sm KK 2 /3+12 s 2 f3 
768 it (3 2 s 2 m KK 2 c w 4 
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Scalar Annihilation into Gauge Bosons: 



Process 


Diagrams 


Total cross section 


zz — ► ZZ 
zz -► W+VF- 

M -> ZZ 
/i/i -> W+iy- 

z0 -> W + W- 

w + w~ -» p^+iy- 

w + w~ — > ZZ 

w + w~ — ► AA 
w + w~ — * AZ 
w + h -► 

-► W+Z 
w+z -> W+Z 


C.3.2 a,b,c,g 
C.3.2 a,b,c,e,f,g 

C.3.2 a,b,c,e,f,g 
C.3.2 a,b,c,e,f,g 

C.3.2 a,b,d,e,f 
C.3.2 a,c,d,e,f,g 

C.3.2 a,b,c,e,f,g 

C.3.2 a,b, c 

C.3.2 a,b,c 

C.3.2 a,b,c,d 

C.3.2 a,b, c,d 
C.3.2 a,b,c,d,e,f,g 
C.3.2 a,b,c,d,e,f,g 


g i (-Asm KK 2 +8m KK 4 )L+if3sm KK 2 +s 2 f3 g 4 H 4, 
64 nc^s' 6 /!' 2 1 512vr s(3c w 4 
(-Wm KK 6 +8m KK 4 s)L-2 /3 s 3 -3 (3 s 2 m KK 2 -S (3 sm KK 4 

64-7T s J (3 2 m,KK 2 
g 4 H 2 (3-2 L gi H 4 
" r 64 C 2,STr/3 2 1 256 I3ksc w 4 
g i (-AsMh 2 +8Mh 4 )L+s 2 f3+4:Mh 2 f3s | l H 4 

64 TTS^C 4 ' 512 /3lTSC w 4 

l (~4(3m KK 2 +8m KK 2 L)H 2 1 8 (3 s+A (3 m KK 2 
512 -k (3 2 sm KK 2 c w 4 1 512 tt (3 2 sm KK 2 c w 4 
g 4 ^ (-4 sMfi 2 +8 M/i 4 )L+s 2 /3+4 M/i 2 /3 s ^ t g i r a 
32 I 7rs a /3 2 J 1 256 /3 7rsc m 4 

9 4 (4m KK 2 c m 2 /3-8m J f X 2 c„ 2 L) J H" 2 g 4 8 c m 4 /3 s+4 c//3 to kk 2 
" l ~256 n /3 2 sm KK 2 c w 4 1 256 i^smffi^ 4 

g 4 5 s 2 /3-48 sm K K 2 i+124 sm if K 2 /9-96 Lm KK 4 
384 tts^ 2 

ff 4 -24m 2 fif sL+24s 2 /3-47smKK 2 /3-4/3m K K 4 
1 768 TTs 2 f3 2 m KK 2 
g 4 17 (3s-96m KK 2 L+l72f3m KK 2 . g 4 H 4 
384 ns 2 /3 2 1 64 /3ttsc w 4 
g i (-3sm KK 2 /3+6sm KK 2 L)H 2 3 s 2 j3+7 sm KK 2 /3-6 sm KK 2 L- f3 m KK 4 


192 irs' 2 f3 2 m KK 2 c w 4 1 7rs 2 /3 2 m K K 2 c» 4 

! g 4 (2c w 2 -i) i (8rn KK 4 L-4m KK 2 Ls+s 2 f3+4sm KK 2 f3) 
64 TTS^C^ 4 

g 4 // 4 S 4 (4)3mM 2 cJ-8m H 2 Lc„, 2 )fi 2 
512s/3ttc w 4 1 512m KK 2 7rs/3 2 c u , 4 
1 8c 1) 4 s/3+4 / 3m mf 2 c™ 4 
~l~ 512m KK 2 nsl3 2 c w 4 
4 (8m JfJf 4 -4sm H 2 )L+4)3sm K 2 4-s 2 /3 

e 27rs 3 /3 2 
9 4 s 2 1 ,(2c 2 -l) 2 (-4sm K 2 +8m K 4 )L+4/3sm K 2 + s 2 /3 

8tt S 3/3 2 

I s w 2 g 4 (l7f3s-96m KK 2 L+l72(3m KK 2 ) 
384 tts 2 /3 2 

X s w 2 g 4 (\7 l3s-9& m KK 2 L+l72f3m KK 2 ) 
384 tts 2 /3 2 

[Eq. (H) below] 
[Eq. (0) below] 



c/ 4 [ -24 /? Ct0 V + 17 s 2 /? c„, 4 + 12 s 2 /? - 96 sc^m^ 2 L + 96 c w 2 sm KK 2 L + 48 sm KK 2 j3 
384 



7T S 3 f3 2 C w 2 

4„„ 2o /IQ„„,_ 2r n«„ 2_ 2o ,no, At „ 2 , nflT^ 4 



172 c w srriKK P — 48 suikk L — 96 c w suikk P — 192 tukk Lc w + 96 Ltjikk 

7T S 3 P 2 C W 2 

60 s 2 /3 c„, 4 — 144 sc w 4 mKK 2 L + 76 c w 4 sm.KK 2 P — 4 c w 4 /? m KK ^ — 48/3 c^ 2 s 2 

27T s 2 f3 2 mKK 2 c w 4: 

-30 c w 2 sm K K 2 l3 + 84 c w 2 sm KK 2 L - 12 sm K K 2 L + 3 sm K K 2 P + 12 s 2 /3 n 



27T s 2 (3 2 mKK 2 c 



2„ 4 



(36) 
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g 4 \-2Af3c w 2 s 2 + 17 s 2 /3c w 4 + 12 s 2 (5 - 96 sc w 4 m KK 2 L + 96 c w 2 sm KK 2 L + 48sm KK 2 /3 
384 [ 7rs 3 /3 2 c „, 2 

172 c w 4 smKK 2 f3 — 48 sm,KK 2 L — 96 c w 2 smKK 2 P — 192 vtikk^LcJ 1 + 96 LniKK 4 

irs 3 (3 2 c w 2 

12s 2 (3- 20 sm KK 2 f3- 4(5 m KK 4 12 H 2 sm KK 2 L - 6 sm KK 2 (3 H 2 3H 4 



2ir s 2 j3 2 mKK 2 



2ir s 2 f3 2 m K K 2 c w 2 



2(3nsc w 4 



(37) 



Scalar Annihilation into Fermions: 



Process 



Diagrams 



Total cross section 



w+z(or(/>) -> /+/- 



C.3.3 
C.3.3 
C.3.3 







2 a 2 
127r/3s 2 



12?r/3s 2 



Here Q and T3 denote the charge and SU(2) charge, respectively, of the fermion. 



D.7 Coannihilation 



D.7.1 Gauge Boson Coannihilation 



Process 



Diagrams 



Total cross section 

(9 z Zq 9imOrn KK '+5s)L-7l0r 



w±wP - 

W±B V - 



(3) 



(3) 



(3) 



-> qq 
* /l/l 

+ /7 



W*£„ M or 
W^B U -► /mu 



C.2.1 a,b 
C.2.1 a,b 
C.2.1 a,b,c 

C.2.1 a,b 
C.2.1 a,b 
C.2.3 a,c,s 
C.2.3 a,c,s,e,d 

C.2.2 a,b,c 
C.2.2 a,b,c 
C.2.2 a,b,c 
C.2.2 a,b,c 

C.2.2 a,b,c,d 



1 



18 



2(AT 2 -1) T2s t tT(P~ 

1 (9^9s 2 )[(10 m ^ 2 + 5 «)^-7 s/3] 

2(iV 2 -l) 72s 2 tt/3 2 
g 4 (5s+14m yK 2 )^+(-16mKX 2 -13s)/3 
576 s 2 7r/3 2 

_1 >/gig 2 ((-10^Ky 2 -55)L+7s/3) 
144 s 2 7r/3 2 

as for -► // 

(-4m KK 2 s-6m KK i )L+2 s 2 (3+5 s/3m KK 2 +ll f3m KK 4 



rriKK' 

g 4 c 2 w {-irn KK 2 s-6m KK 4 ^L+2 s 2 /3+5 s(3 m KK 2 +11 (3 m KK 4 



18 



+ 



injfx 2 ' s 2 /3 2 
g 4 c?u -4m K 2 +s 
576 J^s 1 
g 4 Y h 2 T 3 (h) 2 s w 2 
24tt s/3c w 2 
1 g 4 s w 2 
192 sttB c w 2 

1 g i s w 2 

192 sir B c w 2 
1 ff 4 s™ 2 

192 s7r/3 Cl „ 2 
1 g*(-s+4m KK 2 ) 
"576 7r s 2 S 
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D.7.2 Fermion Coannihilation 



Process 


Diagrams 


Total cross section 


fB v - f(A v + Z v ) 


C.4.1 a,c 


g 2 g 2 Y 2 (T 2 +t a ,ii 2 (9 w )Y 2 )(Y f 2 (6L-8)) 
96tt s/? 2 


JWP - f(A„ + Z„) 


C.4.1 a,c 


g 4 T 2 (T 2 +ta,n 2 (e w )Y 2 ){Ye L 2 (6 L-/3)) 
96tt s/3 2 


J L v J L v 


C 4 1 a r 


9 2 9 2 (F ei 2 (6L-/3)) 
192tt s/3 2 


f~M/ 3 . f+W— 

JL W u ^ Jl W u 


P /I 1 Q K 


c/ 4 (-26 m^if 2 i+32 s/3+23 /3 m K K 2 ) 
768?r s3 2 m KK 2 


W~T, + — > Tr(Z -1-4 "l 


fill q h r 


g 4 ((-32c w 2 +6)m KK 2 L+(24/3c w 2 -8)m KK 2 +32c w 2 sl3) 
768?r s3 2 m KK 2 c w 2 


-> L+(Z„ + A,) 


C.4.1 a,b,c 


as above 


W7<# -QZ(^ + A0 


C.4.1 a,b,c 


g 4 (-6m KK 2 L+240 c ro 2 m ifx 2 L+/3m KK 2 -208/3c tlI 2 m KK 2 -288 Cu , 2 s/3) 


6912c tu 2 7r sf3 2 mKK 2 


C.4.1 a,b,c 


as above 




C.4.1 a,b 


! 9 4 (3m ffi 2 L+4/3s) 
192 TTs8 2 m KK 2 




C.4.1 b,c 


g 4 -16 Lm H 2 +11 f3m KK 2 +8 s/3 
2 1927rs/3 2 m KK 2 


/Z^ + >^/l^ + 


C.4.1 b,c 


as above 




C.4.1 a,c 


9 2 T 3 2 <y 2 6L-/3 
72 sn/3 2 




C.4.1 a,c 


9 2 iYq9 2 s 6L-/3 
72 svr/3 2 




C.4.1 a,c 


9 2 g 2 6L-/3 
36 stt/3 2 


g^q^ qg^ 


C.4.1 a,b,c 


g 4 (25 0m KK 2 -24Lm KK 2 +36 s/3) 
1296vr s3 2 m KK 2 


g^q qAp, z^Wp 


C.4.1 a,c 


9s9%z.w -13+6 L 
2(N 2 -1) 96tts/3 2 
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D.8 Higgs Coannihilation 



JT lUtcoo 








P ^ 1 a h p 


9 4 

96ttc 2 s/3 
o 4 s^ 


hB — ► hZ 


C 5 1 a b c 


y 

96?rc4 s/3 


TT7"3 r7 

zW„ — »■ 2;/„ 


C.5.1 a,b,c 


/ 
967rc^ s/3 




P ^ 1 a h r 


96ttc* s/3 


Mot z)wf> 3) -► w+W- 


C.5.1 a,b,c,d 


g A -4m KK 2 L+/3m KK 2 +4sl3 
96 7r s/3 2 mfcK 2 


0(or -> w+VF,7 


C.5.1 a,b,c 


9 4 sl 
96c^, /3 7r s 


WW, 3 , -»■ + Ajy) 


C.5.1 a,b,c 


<? 4 

967T c w 2 s/3 


^ -► 0W„ 


C.5.1 a,b,c 


4 9 
WD nSfjCw 


™wj 3) -► 0W„ 


C.5.1 a,b,d 


g 4 -4Lm KK 2 +f3m KK 2 +4sP 
96 7r s(3 2 rriKK 2 


<^>(or z)W„ — > 0(or -z)Wj, 


C.5.1 a,b,c 


<? 4 
96?rs/3 


z(or fc)W„ -► /i(or z)W„ 


C.5.1 a,b,d 


j 4 (4s/3+/3m K 2 -4Lm K 2 ) 
96s7r (3 2 rriKK 2 


uTW+->z(or /0(Z„+A„) 


C.5.1 a,b,c,d 


g i (-4Lc w 2 +(3)m KK 2 +4c w 2 sf3 
96 7r s/3 2 c w 2 mKK 2 




C.5.1 b,c,d 


.9 4 


w + W+ -► u> + W+ 


C.5.1 a,c,d 


9 4 (-2L+/3)m M 2 +s(3 
24 7r s(3 2 mKK 2 




C.6.1 a 


(9f94>) 2 sfi-Lm KK 2 

8 7T S0 2 m,KK 2 



In the last line, and denote the couplings of the appropriate ^-channel gauge boson to 
the fermions and the scalars, respectively. 
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